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Abstract. — A systematic review of Euglandina singleyana was undertaken 
to thoroughly examine shell morphology, allozyme, and mtDNA variation in 
specimens of E. singleyana endemic to central Texas. Allozyme similarity 
ranged from 95% in the most geographically proximal individuals of E. sin- 
gleyana (Val Verde and Real counties) to 82% between the most distant indi- 
viduals (Comal and Val Verde counties). DNA sequence similarity, based on 
a 397 bp partial 16S mtDNA sequence, ranged from 98% in eastern specimens 
(Comal and Kerr counties) to 95% in individuals from Kerr, Real, and Val 
Verde counties. Analysis of variation in shell morphology, allozyme similarity, 
and mtDNA sequences supports the existence of a single, highly variable, wide- 
spread species of Euglandina in central Texas. This study also examines the 
validity of E. immemorata, using morphometric and DNA sequence data and 
E. exesa, using morphometric data. The morphometric analysis showed that E. 
immemorata and E. singleyana differ significantly in shape. The current status 
of E. immemorata and E. exesa are also examined. 



Euglandina singleyana (Binney, 1 892) is 
found in a wide variety of habitats along 
the southern margin of the Edwards Plateau 
in Texas, from Terrell County in the west 
to Fayette County in the east, and south to 
Refugio County (Fig. 1; Singley 1893, Pils- 
bry 1946, Fullington & Pratt 1974, Hu- 
bricht 1985). In the eastern part of its range 
it is found under rocks and logs in wooded 
stream valleys in the limestone of the Ed- 
wards Plateau. In the clay and sandy areas 
of the Balcones Escarpment, it is restricted 
to wooded lowlands. The range of this spe- 
cies extends to the western Stockton Pla- 
teau where it is found under fallen Yucca 
and in rock crevices in desert shrub habitat 
dominated by Lechuguilla cactus (Fulling- 
ton & Pratt 1974). 

A great deal of the historical taxonomic 
confusion regarding this species appears to 
be related to its relatively wide geographi- 
cal range of —250 km, compared to the me- 



dian range of land snail distribution of 50 
km reported by Solem (1984) and high lev- 
el of variation in shell morphology. Shells 
collected in the eastern part of the range can 
readily be distinguished from those collect- 
ed in the western part of the range. One 
purpose of this study was to examine vari- 
ation in shell morphology of specimens 
from throughout the range of the species in 
central Texas to determine if there are dis- 
tinct differences in western versus eastern 
shell morphology or if there is continuous 
(clinal) change in shell shape and size 
across the range of the species. The purpose 
of this study was also to examine allozyme 
and DNA sequence variation in individuals 
from across the range of the species (Co- 
mal, Kerr, Real, and Val Verde counties) to 
determine if specimens conforming to the 
description of E. singleyana formed a 
monophyletic group. This study addresses 
the taxonomy of this species from a phy- 
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Fig. 1 . The amended distribution of Euglandina singleyana in Texas. Counties in the historical literature are 
shown in gray. Counties added in this study (Blanco, Chambers, Coryell, Kimble, and Uvalde) are shown in 
black. Collection localities for allozyme and DNA studies are represented by a black square. Museum localities 
listed in the material examined are represented by a black circle. The distribution of this species is not continuous 
throughout the highlighted area, but instead occurs in isolated patches of appropriate habitat. 



logenetic species concept approach using 
both monophyly (sequence analysis) and 
diagnosability (shell morphometries) crite- 
ria (Minton & Lydeard 2003). 

The carnivorous land snail genus Eu- 
glandina has a problematic taxonomic his- 
tory. Von Martens (1901:47) noted that 
“many species have been described only 
from one or a few examples, and not fig- 
ured. In this genus ... it seems to be very 
difficult, or rather impossible, to draw a 
clear line of distinction between local var- 
iations and nearly allied species.” In addi- 
tion to these problems, Thompson (1987) 
also noted that many original descriptions 
are scattered among nineteenth century 



journals in several languages and often lack 
critical details of sculpture of the adult and 
embryonic shells. 

Euglandina singleyana from central Tex- 
as appears typical of this group in having a 
history of confusion regarding its identity. 
Initially, Roemer (1849) identified speci- 
mens from New Braunfels as Glandina 
truncata (Gmelin, 1788). Binney & Bland 
(1869) later considered Texas specimens to 
be Glandina corneola Pfeiffer, 1857. Bin- 
ney (1892) described and named this Texas 
snail Glandina singleyana. Pilsbry later 
(1907) noted that in the Binney (1885) 
monograph, shells from Texas were listed 
with the name G. decussata (Deshayes, 
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Fig. 2. Landmarks used in the geometric morphometric analysis. Numbers 1 and 12 were used as baselines 
in all analysis. Number 1 is the most distal point on the columella to Number 12 at the apex of the shell. 



1850), the description and figure were of G. 
corneola, and the locality and anatomical 
descriptions were those of G. singleyanci. 
Singley (1893) listed this species as G. de- 
cussata var. singleyanci. Pilsbry (1907; p. 
175) removed all Mexican and mainland 
members of Glandina to the genus Euglan- 
dina due to taxonomic confusion regarding 
the name and types of Glandina and sub- 
sequent workers (Pilsbry & Ferriss 1906, 
Pilsbry 1907, 1946, Fullington & Pratt 
1974, Hubricht 1985) have retained this ge- 
neric name. 

Euglandina immemorata Pilsbry, 1907 is 
another species associated with some de- 
gree of taxonomic uncertainty. This species 
was described by Pilsbry (1907) on the ba- 
sis of two shells from “Texas” (exact lo- 
cality unknown). Fullington & Pratt (1974) 
stated that the holotype of E. immemorata 
is not especially distinct from E. singleyanci 
but could not be precisely matched by any 
material they had seen. They concluded that 
E. immemorata is probably E. singleyana 
or, if distinct (based on the improved 
knowledge of the Texas fauna since 1946) 
is probably not a Texas species. 

Euglandina exesa Cockerell, 1930 was 
described from a single shell found in a 
limestone deposit in a cinnabar mine at Ter- 
lingua, Brewster County, Texas. Fullington 
& Pratt (1974) state that the holotype does 



not differ from many western specimens of 
E. singleyana. They consider E. exesa to be 
simply a western range extension of E. sin- 
gleyana. 

A systematic review of E. singleyana 
was undertaken to thoroughly examine the 
observed morphological differences and 
add further evidence using genetic tech- 
niques. This study also tests the conclusions 
of Fullington & Pratt (1974) by examining 
the validity of E. immemorata using mor- 
phometric and DNA sequence data and E. 
exesa using morphometric data. 

Materials and Methods 

Morphometric analysis . — Programs used 
for morphometric analysis are part of the 
Integrated Morphometries Package (www. 
canisius.edu/~sheets/morphsoft.html) made 
available by David Sheets, Miriam Zeld- 
itch, and Donald Swiderski. Specimens of 
Euglandina singleyana (64), E. immemor- 
ata (six — including both type specimens), 
and E. exesa (holotype) were examined 
(Appendix 1). Twelve landmarks on each 
shell were digitized from photographs using 
the program tpsDig version 1.31 (F. J. 
Rohlf. tpsDig32: Digitize coordinates of 
landmarks and capture outlines. http.V/life. 
bio.sunysb.edu/moiph/index.html) (Fig. 2). 
The landmarks were chosen to be repeat- 
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able and homologous (Swiderski 1993) 
across all shells, and most landmarks rep- 
resent shell sutures (Stone 1998). For each 
specimen, digitized landmark coordinates 
were transformed to procrustes distances 
using partial procrustes superimposition 
methods. This was carried out in Coord- 
Gen6 (H. D. Sheets, Dept, of Physics, Can- 
isius College, 2001 Main St. Buffalo, New 
York 14208, sheets@gort.canisius.edu). A 
principal component analysis was per- 
formed using partial procrustes distances 
(PCAGen6e; H. D. Sheets). This program 
computes partial warp scores for each spec- 
imen, using a procrustes mean specimen 
based on all data for comparison. The prin- 
cipal components (eigenvectors of the co- 
variance matrix) are then calculated based 
on the covariance matrix derived from the 
partial warp scores. Landmarks 1 and 12 
were designated end points and all speci- 
mens were standardized according to this 
baseline. This transformation to partial 
warp scores accomplishes standardization 
so that the principal component analysis ex- 
amines variation in shape, excluding vari- 
ation due to scale, rotation, and translation 
(Swiderski 1993, Stone 1998). 

A canonical variates analysis was per- 
formed using the program CVAGen6 (H. D. 
Sheets) to determine the set of axes that al- 
lows for the greatest possible ability to dis- 
criminate between two or more groups. 
This program computes partial warp scores 
with reference to a common mean specimen 
then performs a multivariate analysis of 
variance followed by a canonical variates 
analysis. It determines how many distinct 
axes there are in the data, (p = 0.05) and 
computes the canonical variates scores of 
all the specimens. It also uses Mahalanobis 
distances to assign specimens to one of the 
groups. The single specimen of E. exesa 
was coded as an unknown and assigned to 
a group based on the canonical variates 
axes formed in the prior analysis (CVA- 
Gen6; H. D. Sheets). Lastly, a comparison 
was performed to determine if there was a 
significant difference in shape between E. 



singleyana and E. immemorata (Two- 
Group6c; H. D. Sheets). 

Allozyme analysis . — Each collection lo- 
cality of Euglandina singleyana was rep- 
resented by one or two specimens. A total 
of seven specimens from five collection lo- 
calities at —50 km intervals throughout the 
geographic range of this species were ex- 
amined (Appendix 1). 

After collection, individuals were held 
without feeding for 7—10 days, then frozen 
in cryotubes in liquid nitrogen and stored 
in an ultracold freezer (— 80°C) until anal- 
ysis. Samples were homogenized in two 
volumes of distilled water using a glass rod 
and centrifuged to obtain an aqueous ex- 
tract. Procedures for cellulose acetate elec- 
trophoresis and staining followed those of 
Hebert & Beaton (1993). Gels were pur- 
chased from Helena Laboratories Inc. 
(Beaumont, Texas), and the buffer used was 
tris-glycine pH 8.5. To examine variation 
within E. singleyana, scorable data for 19 
loci (Table 1) were obtained and analyzed 
using Tools for Population Genetic Analy- 
sis 1.3 (M. P. Miller. Tools for Population 
Genetic Analysis (TFPGA 1.3): A windows 
program for the analysis of allozyme and 
molecular population genetic data. Com- 
puter software distributed by the author.). 
To determine genetic similarity, Nei’s un- 
biased genetic identity was calculated (Nei 
1978). An unweighted pair group method 
using arithmetic averages (UPGMA) cluster 
analysis was then performed using the ge- 
netic identity matrix. 

DNA sequence analysis . — Twelve tissue 
samples were either preserved in 70% eth- 
anol or frozen in liquid nitrogen and then 
stored in an ultracold freezer at — 80°C. The 
outgroup taxa chosen were the closest rel- 
atives of Euglandina with sequences avail- 
able on GenBank, relationships from Wade 
et al. (2001). Total genomic DNA was ex- 
tracted from several milligrams of foot tis- 
sue by digestion with lysis buffer and Pro- 
teinase K and then purified by phenol: chlo- 
roform extraction according to standard 
procedures (see Palumbi, S., A. Martin, S. 
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Romano, W. O. McMillan, L. Stice, & G. 
Grabowski. 1991. The Simple Fool’s Guide 
to PCR. Privately distributed, Honolulu, 
Hawaii, 40 pp.). 

Mitochondrial DNA sequences were ob- 
tained for an amplified segment of the 16S 
rDNA gene using 16sar and 16sbr primers 
(see Palumbi above). Approximately 10 ng 
of genomic DNA provided templates for 
double-stranded reactions via the polymer- 
ase chain reaction (PCR). PCR reactions 
were done in a 50 |xL solution containing 
each dNTP at 0.22 pM, each primer at 0. 1 
|xM, 1.5 mM MgCl 2 . 1 unit Taq DNA poly- 
merase, and IX PCR reaction buffer. Re- 
actions were amplified for 30 cycles of 
92°C for 45 sec, 50°C for 45 sec, and 68°C 
for 2 min. Samples were purified and dou- 
ble-stranded DNA provided the template 
for cycle-sequencing using BigDye (ABI) 
chemistry followed by analysis on an 
ABI3100 automated sequencer. 

Contigs were assembled in Sequencher® 
4.0.5 (Gene Codes Corporation, Ann Arbor, 
Michigan and aligned by eye using BioEdit 
(Hall 1999) with reference to secondary 
structure models to refine the alignment and 
identify regions corresponding to loops and 
stems (Lydeard et al. 2000). Sequences were 
deposited with Genbank (Accession Num- 
bers: AF405235-AF40524 1 , AY149279, 
AY167887— AY167889). Aligned sequences 
were analyzed using maximum parsimony 
with PAUP*4.0bl0 (S wofford 2002) using a 
heuristic search (10 addition replicates). The 
following options were used: uninformative 
characters were ignored, only minimal trees 
were kept, gaps were treated as missing, and 
zero length branches were collapsed. A 
bootstrap analysis with 1000 iterations was 
conducted (Felsenstein 1985). Bremer sup- 
port values (Bremer 1994) were calculated 
using the Decay function of MacClade 4.03 
(Maddison & Maddison 2000). 

Results and Discussion 

Morphometric analysis. — In the exami- 
nation of shell variation within Euglandina 



Table 1. — Presumptive enzymatic loci resolved. 
Tris-glycine (pH 8.5) was used as the buffer system 
for all of the enzymes listed. All stain and buffer rec- 
ipes are from Hebert & Beaton (1993). 



Enzyme system and abbreviation 


Enzyme 

commission 

number 


Adenylate Kinase (ADK) 


2.1 A3 


Alcohol Dehydrogenase (ADH) 


1.1.99.8 


Alcohol Oxidase (AOX) 


1.1.3.13 


Aspartate Aminotransferase (AAT) 


2.6.1. 1 


Glucokinase (GK) 

Glucose-6-phosphate Dehydrogenase 


2.1 A. 2 


(G6PDH) 


1.1.1.49 


Hexokinase (HK) 


2.7.1. 1 


L-Iditol Dehydrogenase (IDDH) 


1.1.1.14 


Isocitrate Dehydrogenase (IDH) 


1.1.1.42 


L-lactate Dehydrogenase (LDH) 


1.1.1.27 


Malate Dehydrogenase (MDH- 1 & 2) 
Malate Dehydrogenase (NADP+) 


1.1.1.37 


(MDHP-1 & 2) 


1.1.1.40 


Nucleoside Phosphorylase (NSP) 


2.3.2. 1 


Phosphoglucomutase (PGM) 
Phosphogluconate Dehydrogenase 


5.4. 2. 2 


(PGDH) 


1.1.1.44 


Superoxide Dismutase (SOD) 


1.15.1.1 


Triosephosphate Isomerase (TPI) 


5.3.1. 1 



singleyana the first principal component 
(PCI) accounted for 67.42% of the varia- 
tion present in the measurements and PC2 
accounted for 9.66% of the variation pres- 
ent. A scatterplot, with specimens grouped 
by county, comparing the first two principal 
component axes, does not show any distinct 
groups within E. singleyana (Fig. 3). Shell 
variation described by principal component 
1 appears to be continuous and the results 
of this analysis do not allow E. singleyana 
to be separated into groups based on shell 
characteristics. 

The canonical variates analysis (Fig. 4) 
to determine if E. immemorata could be 
distinguished from E. singleyana was sig- 
nificant, (df = 20, p = 3.70685 X 10~ 5 ) 
with 100% correct reclassification of both 
species. A two-group comparison showed 
that the shapes of the two species are sig- 
nificantly different (Hotelling’s T 2 /= 6.52, 
df= 24, 45, p = 3.8316 X HF 8 ). This result 
indicates that E. immemorata is morpholog- 
ically distinct from E. singleyana. This re- 
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Fig. 3. Scatterplot displaying principal component scores of Euglandina singleyana. X-axis is principal 
component one, Y-axis is principal component two. Texas counties are represented by the following symbols: 
Bexar (asterisk) type-locality. Chambers (circle), Comal (plus sign), Coryell & Travis (star), Hays (square), 
Kimble, Uvalde, & Real (hollow triangle), Val Verde (filled triangle), Victoria (X). 



suit, along with the discovery of specimens 
of E. immemorata in Nuevo Leon, Mexico 
(Correo-Sandoval 1993), supports the re- 
moval of E. immemorata from Texas faunal 
listings. 

The unknown assignment test to deter- 
mine the placement of E. exesa placed the 
single specimen in a cluster formed by E. 
singleyana (Fig. 4); however, this result was 
not significant (p > 0.001). This species is 
represented by only one specimen. There- 
fore, this test is not replicable and has little 
statistical power. It is interesting that the 
CVA placed the specimen of E. exesa with- 
in E. singleyana, supporting the statements 
of Fullington & Pratt (1974) about this spe- 
cies. However, without more evidence, 
nothing conclusive can be stated about the 
validity of this species. 

Allozyme analysis . — The results of this 
study reveal a moderate to high degree of 



genetic similarity in all specimens of Eu- 
glandina singleyana examined across its 
range. Seven of 19 loci (MDH-2, MDHP- 
2, IDH, LDH, SOD, G6PDH, PGM) were 
monomorphic for all specimens. The Val 
Verde County specimens exhibited poly- 
morphisms at nine enzymatic loci (MDH-1, 
MDHP-1, PGD, AOX, A AT- 1 , ADK, 
IDDH, ADH, HK), the Real County speci- 
mens exhibited polymorphism at six 
(MDH-1, MDHP-1, AOX, HK, NSP, 
ADH). In the specimens examined, there 
were no fixed allelic differences among 
populations of E. singleyana across its geo- 
graphic range. 

Allozyme variation appeared to follow a 
general geographic pattern with specimens 
having the smallest genetic distance from 
their most proximal geographic neighbors 
(Fig. 5). Genetic divergence calculated us- 
ing Nei’s unbiased genetic identity resulted 
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Fig. 4. Graph of the canonical variates axes displaying the separation between the means of Euglandina 
singleyana and E. immemorata. Plotted as an unknown is E. exesa . Oversized symbols indicate the means. X- 
axis is CVA 1, Y-axis is CVA 2. 



in the most western specimens from Val 
Verde County and the Real County speci- 
mens displaying a genetic similarity of 
94.5%. The Val Verde + Real cluster was ge- 
netically similar to the more centrally lo- 
cated Kerr County specimen at a level of 
89.7%. This cluster is 85.2% similar to the 
most eastern specimen from Comal County. 

Perez & Strenth (2002) found that spec- 
imens of Euglandina texasiana (Pfeiffer, 
1 857) from collection localities 150 km dis- 
tant, in south Texas and northern Tamauli- 
pas, had a similarity of 94.5%. In compar- 
ison, this study observed 94.5% similarity 
between specimens of E. singleyana located 
157 km apart, and 85.2% similarity be- 
tween the most geographically distant spec- 
imens (246 km) from Comal and Val Verde 
Counties. Euglandina singleyana was 
found to display enzyme polymorphism at 



12 of 19 loci examined. This level of allo- 
zyme variability is more similar to the re- 
sults of studies on Helix aspersa (Selander 
& Kauffman 1975) and differs from results 
found in Liguus by Hillis et al. (1991) and 
Rumina decollata by Selander & Kaufman 
(1975), which were notably less variable. 

The levels of genetic distance among the 
widely separated populations of Euglandina 
singleyana are within the range found in 
other organisms for genetic divergence be- 
tween subspecies (Quicke 1993). However, 
there are no fixed differences among the 
populations and genetic distance is low be- 
tween geographically proximal populations. 
This analysis indicates that the specimens 
of E. singleyana examined represent a sin- 
gle species. 

DNA sequence analysis . — The region of 
16S mtDNA that was sequenced resulted in 
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Fig. 5. A dendrogram of allozyme similarity in Euglandina singleyana produced by UPGMA on a matrix 
produced using Nei’s unbiased genetic identity. 



an aligned data matrix of 397 base pairs of 
which 191 were phylogenetically informa- 
tive. Phylogenetic analysis of the data using 
maximum-parsimony analysis resulted in 
three equally parsimonious trees of 480 
steps (Cl = 0.8646, gl = -2.024336). 
These three trees differed only in the place- 
ment of the Comal County (eastern) speci- 
men relative to the Kerr County specimen. 
One topology places the Kerr County spec- 
imen (central) sister to the Real + Devil’s 
River clade (Central & Western). An alter- 
nate topology places the Comal County 
specimen in this position, and the final to- 
pology describes Comal and Kerr County 
specimens as each others’ closest relative. 
Bootstrap analysis (1000 pseudoreplicates) 
of the aligned data matrix using maximum- 
parsimony produced the tree shown in Fig. 
6. Euglandina singleyana forms a mono- 
phyletic group. The sister taxon to E. sin- 
gleyana is the specimen of Euglandina 
from Northern Coahuila. Also outside this 
grouping is E. corneola from Tamaulipas, 
Mexico. Pairwise sequence identity was 
calculated for each clade. Within the Val 
Verde County cluster (Devil’s River and 
Comstock) there is 98—100% sequence sim- 
ilarity. The two Real County specimens had 
identical 16S sequences (=100% similari- 
ty). Sequence similarity between Ken- 



County and Comal County was 97%. All 
possible combinations of pairwise compar- 
isons were performed with a minimum sim- 
ilarity of 95% among specimens from Ken 
and Real County. The two specimens of E. 
immemorata formed a group apparently not 
closely related to E. singleyana. 

Both allozyme and sequence analyses 
show an interesting geographic pattern with 
populations most closely related to their 
geographically proximal neighbors. This 
pattern of strong geographic structuring is 
often seen in land snail species (Thomaz et 
al. 1996, Schilthuizen et al. 1999). Thomaz 
et al. (1996) examined geographic variation 
within Cepea nemoralis and Helix aspersa 
and found very high levels of sequence di- 
vergence (12%) within these species of land 
snails. The authors conclude that the most 
likely explanation for the observed levels of 
divergence is the population structure of 
land snails with low dispersal and large 
populations divided into infrequently inter- 
acting demes. 

Avise et al. (1987) presents several tests 
of this hypothesis, one of which is: phylo- 
genetic differentiation between long sepa- 
rated demes should be reflected in nuclear 
as well as mitochondrial assays. The con- 
gruence of the allozyme and mitochondrial 
data in the present study appear to fulfill 
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Fig. 6. Strict consensus of three most parsimonious trees of partial 16S mtDNA sequences of specimens of 
Euglandina singleyana from Comal, Kerr, Real, and Val Verde counties. Numbers above the branches are 
bootstrap support percentages (1000 pseudoreplicates). Numbers below the branches are Bremer’s support values. 



this requirement, lending credence to the 
idea that the specimens of E. singleyana ex- 
amined in this study represent a single spe- 
cies. The present study does not have 
enough samples to conclusively examine 
geographic partitioning within this species; 
however, there does appear to be some geo- 
graphic structure in the data. 

Although sample sizes were low, they 



appear to be adequate to address the ques- 
tions posed in this research. Euglandina 
singleyana is considered to be an uncom- 
mon species (Singley 1893; Neck 1984, 
1988), and the area where they are found is 
rapidly being disturbed ecologically due to 
a growing human population. As a result of 
their very specialized feeding habits, rarity, 
habitat preferences, and human activities. 
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the live specimens examined during this 
study represent a significant collection of 
living specimens of Euglandina singleyana. 

Conclusions 

All of the analyses, including morpho- 
metric, allozyme, and mtDNA sequences, 
support the premise that Euglandina sin- 
gleyana is a single, widespread, highly var- 
iable species. Both allozyme and DNA se- 
quence results indicate that there are some 
detectable genetic geographic patterns with- 
in E. singleyana, as well as the observable 
morphological gradient across the range of 
this species. The geographically most dis- 
tant specimens are the most genetically dis- 
similar and geographically proximal speci- 
mens are more similar. 

The results of the morphometric and se- 
quence analyses indicate that E. immemor- 
ata is distinct from E. singleyana. This re- 
sult, combined with the recent discovery of 
specimens of E. immemorata from Nuevo 
Leon, Mexico leads to the conclusion that 
E. immemorata should be removed from 
Texas faunal listings. The taxonomic place- 
ment of Euglandina exesa could not be de- 
finitively determined. 
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Appendix 1 

Material Examined and Distribution 
Morphometric Analysis 

Euglandina singleyana (64): The number of speci- 
mens of each lot that were used follows the museum 
catalog number in parenthesis (all undamaged shells 
with 5 or more whorls were used). Field Museum of 
Natural History: 29 specimens: FMNH 22347 (1), 
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30850 (2), 36309 (5), 5001 1 (4), 58295 (3), 62416 (6), 
74849 (1), 78560 (2), 98199 (2), 109097 (1), and 
175899 (1). The Academy of Natural Sciences of Phil- 
adelphia: 22 specimens: ANSP 186729 (1), 186731 
(1), 346372 (3), 84622 (2), 104753 (1), 4312 (1), 
84637 (1), 158379 (1), 134180 (4), 76849 (3), 150798 
(1), 76837 (1) topotype, 186733 (1), and 87425 (1). 
Angelo State Natural History Collection: 13 speci- 
mens, individually numbered: (ASNHC 0008-0014, 
042, 044-048). 

Euglandina immemorata (6): FMNH 11777 (1), 
4359 (1) locality unknown, types; Florida Museum of 
Natural History 189621 (1), Nuevo Leon, Mexico, 
Santiago, 1 km. N of Laguna de Sanchez. University 
of Alabama Gastropod Collection, 3 specimens, indi- 
vidually numbered: 632-634, Nuevo Leon, Mexico 
25°23'00.9"N; 100°14'28.6"W, 1 km N of Laguna de 
Sanchez. 

Allozyme Analysis 

Euglandina singleyana (7): Specimens used in this 
analysis are deposited in the Strecker Museum of Nat- 
ural History, (SMNH) Baylor University (Accession 
Number 2001-A-l-l; Catalog Numbers SM32439- 
SM32446). (1) Landa Park, New Braunfels, Comal 
County, Texas; (1) Kerrville, Kerr County, Texas; (2) 
9 miles N of Leakey, Real County, Texas; (2) Devil’s 
River State Natural Area (DRSNA), Val Verde County, 
Texas; (1) Comstock, Val Verde County, Texas. Rab- 
dotus alternatus from Val Verde County was selected 
as an outgroup. 

DNA Sequence Analysis 

The same specimens listed in the above section on 
allozyme analysis provided the tissue samples which 



were used in the sequence analysis. Additionally in- 
cluded were a specimen of an unidentified Euglandina 
from La Cuesta in Northern Coahuila; Euglandina cor- 
neola from 4 miles SW of Mante, Tamaulipas, Mexico; 
2 specimens of E. immemorata from Nuevo Leon, 
Mexico 25°23'00.9"N; 100°14'28.6"W, 1 km N of La- 
guna de Sanchez; and one specimen of E. rosea from 
Lake County, Florida 28°30'53"N; 81°44'15"W. 

Euhadra amaliae (AF098712), Phaeohelix phaeo- 
gramma (AF098714), and Nesiohelix omphalina 
(AF098713) sequences from Chiba (1999) were used 
as outgroups. 

Distribution 

Strecker (1935) listed Euglandina singleyana as oc- 
curring in Bexar, Caldwell, Comal, Goliad, Hays, 
Travis, Victoria, Atascosa, Frio, Gonzales, Guadalupe, 
and Wilson counties. McGee (1971) added Bandera, 
Fayette, Kendall, Kerr, Kinney, Val Verde, Medina, 
and Real counties. Cheatum et al. (1972) added Terrell 
County to the distribution. Collections by the authors 
and examinations of museum specimens from the 
Academy of Natural Sciences in Philadelphia and the 
Field Museum of Natural History during this study add 
Uvalde, Kimble, Chambers, Coryell, and Blanco coun- 
ties to the known distribution (Fig. 1). 

County records. — FMNH 78560, Horse Creek, Cor- 
yell Co., Texas; FMNH 98199, Galveston, Chambers 
Co., Texas; FMNH 210, 103/1, River Drift, Pedernales 
Falls State Park, Blanco Co., Texas; ANSP 186731, 
Garner State Park, Uvalde Co., Texas; ANSP 186727, 
Roadside Park, Nueces River W of Uvalde, Uvalde 
Co., Texas; ANSP 186733, Llano River, Highway 29, 
12 miles S of Junction, Kimble Co., Texas. 



